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We develop the hydrodynamic theory of Fermi superfluids in the presence of a periodic potential. 
The relevant parameters governing the propagation of sound (compressibility and effective mass) are 
calculated in the weakly interacting BCS limit. The conditions of stability of the superfluid motion 
with respect to creation of elementary excitations are discussed. We also evaluate the frequency of 
the center of mass oscillation when the superfluid gas is additionally confined by a harmonic trap. 
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It is well known that sound cannot propagate in a de- 
generate noninteracting Fermi gas. In fact the sound 
velocity c = vp/\/3, calculated from the compressibil- 
ity of the gas, is smaller than the Fermi velocity vp, so 
that sound is unstable towards decay into the continuum 
of particle-hole excitations. In the presence of an attrac- 
tive interaction, Cooper pairing opens a gap at the Fermi 
level. As a result, sound waves with sufficiently low fre- 
quency can propagate in a neutral superfluid. An impor- 
tant question is to understand how these modes behave 
when the system is confined by a periodic potential. 

The propagation of sound in Bose gases trapped by an 
optical lattice has been already the object of extensive 
theoretical work 1]. Experiments on the collective 
oscillations have confirmed that the low energy dynamic 
behaviour of a Bose-Einstein condensed gas is well de- 
scribed by the hydrodynamic equations of superfluids ^ , 
where the coherent tunneling through the barriers gen- 
erated by the periodic potential is taken into account by 
a proper renormalization of the effective mass. 

In this paper we generalize the hydrodynamic equa- 
tions of Fermi superfluids in order to investigate the prop- 
agation of sound and the collective oscillations of a Fermi 
gas in the presence of a ID periodic potential generated 
by a laser field. 

We consider an interacting two-component atomic 
Fermi gas trapped by a potential Vext given by the sum 
of a 3D harmonic trap and of a stationary ID optical 
potential modulated along the z-axis: 



Vho + sEr sin qbz, 



(1) 



where Vho = m {uj±r'j_ + uj^z"^) /2 and Er = h'^q'^/2m is 
the recoil energy. Here qb is the Bragg momentum and 
s is a dimensionless parameter providing the intensity 
of the laser beam. The optical potential has periodicity 
d = tt/qb along the z-axis. 

The long wavelength dynamic behaviour of a neutral 
superfluid is properly described by Landau's hydrody- 
namic theory. At T = the relevant variables are the 



density n of particles and the superfluid velocity 



(2) 



where m is the atomic mass and (j) is the phase of the 
order parameter defined through the anomalous mean 
value (vI/T(r)vI/x(r)) = |(vI'^(r)M/4(r))|e'*«. 

The hydrodynamic equations can be derived from a 
variational principle, starting from an effective action 
describing the low energy collective excitations over the 
ground state. In the weak coupling limit these excita- 
tions are known as the Bogolubov- Anderson modes. The 
Lagrangian of the system takes the general form 



L= dr 



e(7i,v) 



(3) 



where e(n, v) is the energy per unit of volume. By ex- 
panding e(n, v) up to the terms quadratic in v one finds 



e(n, v) = e(n, 0) + ]:nw\ + ]- ^nv^ (4) 
2 Z m 

where rh (> m) is an effective mass accounting for the 
increased inertia of the superfluid along the direction of 
the laser. In general rh is density dependent. In Eq.Q 
we have assumed that the two spin species are equally 
populated (n = 2n|) . Inserting Eq.Q) in the Lagrangian 
(PI yields the hydrodynamic equations 

— + dx{nvx) + dy{nvy) + dz{—nvz) = 0, (5) 



m m o d ,m ^ m ^ 
'-Vho+^vl + ^vl + — (-n)-vl 



0, 



— —\7( 

(6) 

where — de/dn is the chemical potential. General- 
ization to include 2D or 3D optical lattices is trivial. 

Let us first consider the propagation of sound in the 
absence of the harmonic potential [Vho = 0). We will 
include later the harmonic trap to investigate the oscil- 
lations of the cloud. If Vho = Eqs. (0) and © admit 



2 



linearized solutions around the equilibrium density no of 
the form n = Uq + Sn, with 6n ^ ^^qz-iujt ^^^^ analo- 
gously for Vz- This permits to derive the phonon disper- 
sion uj = CzQ along the direction of the optical lattice, 
with the sound velocity given by 



n dfj, 
ifi dn 



(7) 



The presence of the periodic lattice enters Eq. (|7|) through 
both the renormalization of rh and the change of the 
inverse compressibility ndfi/dn. Analogously, the sound 
velocity along the transverse directions, where the gas is 
free, is given by — {n/m)d^/dn. 

Equations (|5I7|I are very general and apply also to 
strongly interacting superfluids. They permit to calcu- 
late the low energy dynamics of the system once the equa- 
tion of state and the effective mass are known. In the fol- 
lowing we will apply these equations to a dilute Fermi gas 
interacting with negative scattering length (BCS limit). 
In this case the internal energy is fixed by the quantum 
pressure and therefore the chemical potential is equal 
to the Fermi energy {fi — ep)- For simplicity we con- 
sider values of the Fermi energy such that only the low- 
est Bloch band is populated at zero temperature. The 
single particle energy spectrum can then be written as 
e(k) ~ k^/2m + e(fcz) where e{kz) is the dispersion re- 
lation of the lowest Bloch band. The energy of the sys- 
tem at rest is given hy E — V J e (k) rii^dk/ {27Th)^, where 
nic = 20 (e^^ — e(k)) is the quasi- momentum distribution 
of the non- interacting gas 4], Q(x) being the usual step 
function, V is the volume of the system and the inte- 
gration over kz is restricted to the first Brillouin zone 
—nh/d < kz < Trh/d. The Fermi energy ep is related 
to the atomic density by the normalization condition 
n = / nkdk/{2Trh)^. 

Let us first assume that the gas moves with constant 
superfluid velocity v. The order parameter acquires the 
phase (j) = 2mv • r/h (see Eq.®) and, correspondingly, 
the field operator ^„ transforms according to 



r/h 



(8) 



Equation (jS)) is a gauge transformation which does not 
change the energy levels, but changes their classifica- 
tion. In particular the Bloch state of quasi-momcntum 
k is mapped into the state k + mv (and hence e(k) into 
e(k TOv)). The energy density of the moving system 
can therefore be written as: 



e(n, v)= / e (k -I- TOv) nk 



(9) 



where the density dependence comes from the quasi- 
momentum distribution nk. We can now expand Eq.Q 
in powers of v. The linear term vanishes by symmetry 
as required by the stability of the ground state. The first 
non vanishing contribution is the quadratic, diamagnetic 



term. Comparison with Eq.Q) then permits to find, after 
an integration by parts, the useful expression 



1_ ^ S{de/dkzf 5{tF 
m 



£(k))dk/(27r?i)3 



n/2 



(10) 



for the effective mass, holding for a dilute superfluid 
Fermi gas at zero temperature. From Eqs Q and 
(|10|l . we also see that the current density along z in 
the equation of continuity (j^l corresponds to jz = 
(l/mV)5E/5vz. 

From Ea. H10() we see that only the states near the 
Fermi surface contribute to the effective mass. In gen- 
eral rn depends on the Fermi energy or, equivalently, 
on the density. An important exception occurs when 
the Bloch band is almost empty. In this case the oc- 
cupied states obey an approximate quadratic dispersion 
e{kz) — fc^/2ni*, where m* depends only on the laser 
intensity. From Ea. HlOjl we then find m = m*. This co- 
incides with the result for the effective mass holding for 
a dilute Bose-Einstein condensed gas in the same optical 
lattice j^l. Equation (|1U|I also shows that the integra- 
tion along the transverse directions plays a crucial role 
in keeping rh finite as the Fermi energy crosses the Bloch 
bandwith. In fact if we neglect the radial dispersion and 
consider a pure ID system, the effective mass (|1(JII would 
actually diverge as the band is completely filled. 

A second crucial ingredient needed to calculate the 
velocity of sound (see Eq.|0) is the compressibility 
{nd^/dn)^^ . In a dilute Fermi gas at zero temperature 
this is simply calculated in terms of the density of states 
as 



j6{eF-e{k))dk/{2nhf 



(11) 



The optical lattice modifies the compressibility with re- 
spect to the value (37r^)^/'^/i^n^/^/3TO holding for a uni- 
form gas. By inserting Eqs (|10l) and Hll|l into Eq.Q, we 
find the result 



J{de/dk,)^5{eF - e(k))dk 
J 6{eF ~ e{k))dk 



(12) 



for , which can be regarded as the average of the square 
of the group velocity de/dkz over the Fermi surface. 
Again it should be noted that, due to the free disper- 
sion in the transverse direction, the sound velocity H12|l 
remains finite as crosses the bandwidth. In a pure 
ID system, Ea. (jl2(l would predict Cz = when the band 
is completely filled, since de/dk^ vanishes at the edge of 
the Brillouin zone. 

Let us point out that Ea. H12l) has been derived start- 
ing from the hydrodynamic equations but the same re- 
sult can be obtained by generahzing BCS theory through 
the inclusion of the interaction between Bogolubov quasi- 
particles. This can be achieved, for example, using time 
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dependent Hartree-Fock theory ^ or the Random Phase 
Approximation (RPA) |7|. 

If the laser intensity is sufficiently large one can work 
in the tight-binding approximation where the dispersion 
of the lowest Bloch band takes the simple form e(fcz) = 
S(l — cos{kzd/h)), with S proportional to the tunneling 
rate between two consecutive wells. In this case the equa- 
tion of state can be calculated analytically. For values of 
the Fermi energy above the bandwidth (e^ > 2(5) one 
finds n = m{eF — 5)/'Kh^d and hence dn/d^ = m/irh^d. 
For smaller values {ep < 2(5) one instead finds n = 
(2^/y(T^^^+{2y-l) ai ccosil - 2y)) m5 / ir'^h^ d, yielding 
dn/dfj, = {m/ir^h^d) arccos(l — 2y), where y = eF/25. 

The effective mass (jlOII increases as the density of par- 
ticles increases. For values of < 25 one finds ffi = 
{2TT'^h'^nl5'^dm)[2{2y ~ l)y'y{\ -y) + arccos(l - 2y)]-^ 
and in the limit y <C 1, corresponding to ei? <C 2(5, m 
approaches the density independent value m* = /5d^. 
Conversely, for tp > 25 the effective mass scales linearly 
with the density: m = 2h'^im/m5'^d. This result can be 
obtained in a BCS approach, following Ref.|3|- 

Let us now use Eq. Q to evaluate the behaviour of the 
sound velocity along the direction of the optical confine- 
ment. In the limit ep > 25 one finds the density indepen- 
dent value = 5^d^ /2'h?. This is simply understood by 
noticing that de/dkz — {5d/h) sm{kzd/h) and that the 
average of sin^^kzd/h) over the first Brillouin zone gives 
the factor 1/2. For ep < 25, the longitudinal sound ve- 
locity is related to the Fermi energy by the expression 



,5d^2l{2y-l)yW 



arccos(l — 2y) 



(13) 



In particular, for ?; ^ 1, corresponding to ei? <C 25, 
Eg. lfT^ reduces to = 2eF/3m*. It is important to 
notice that the continuous hydrodynamic approach along 
the direction of the laser field requires that the phase cj) 
of the order parameter varies slowly between neighboring 
wells. Taking into account that the energy of the phonons 
must be small compared to the superfluid gap A, we see 
that our theory is valid provided ^ m.m{A/ he z, d~^) 
where Qz is the wave vector of the sound excitation. 

The propagation of sound in the transverse direction 
is affected by the presence of the optical lattice only 
through the change of the compressibility. For values 
of the density corresponding to > 25, one finds 
c\ = (ep — 5)/m. For smaller values of one finds 
c\ = TT^ dn / a.TCCOs{l — 2y), yielding c\ — 2ep/3m 
in the limit ep <^ 25. In the case of transverse sound 
propagation, the hydrodynamic theory is valid under the 
usual assumptions q± <C A/hc±. 

The hydrodynamic equations (|5I6I) can be generalized 
to include finite velocity of the gas with respect to the 
optical lattice. This is important in order to study the 
conditions of stability of the moving fiuid with respect to 
the creation of elementary excitations. This problem has 



been recently addressed in the case of Bose-Einstein con- 
densates, both from the theoretical and experimental 
01 point of view. 

Starting from the general form Q for the Lagrangian 
and setting Vho = 0, one finds 



— +Vj(n,v) =0, 
m— + V^(n,v) = 0, 



(14) 
(15) 



where j = {l/m)de/dv, /j, = de/dn are, respectively, 
the current and the chemical potential of the gas while 
e(n,v) is defined in Eq.jnj. 

Linearizing Eq.s (|14I15|I with n(r) — no + 5n{z) and 
v(r) — [k/m + 5v{z))i and looking for solutions of the 
form 5n, 5v ^ e*^*^~"*\ we find the dispersion relation 



dndk ' 



I9l 



(16) 



Dynamical instability occurs when the argument in the 
square root of Ea. (|16|l becomes negative. This leads to 
complex frequency and corresponds to an exponential 
growth of excitations wich destabilize the superfluid. En- 
ergetic instability instead takes place when lo becomes 
negative, i.e. when the first term in Ea. H15|l (account- 
ing for the Doppler shift) is in modulus larger than the 
second one. 

In tight-binding limit the energy density takes the 
simple form 



t(n, k) — e{n, 0) 



(17, 



where m is defined in Ea. (|10|l . Note that in the limit 
of smaU k Ea. lfTTjl reduces to Eq.Q with = 0. By 
inserting Ea. (|17(l into Ea. H16|l we find 



d / n\h 
dn \rhJ d 



kd 



d'^e n rkd^, , 

— -^COs[—]\q\. (18) 

on^ m n 



We see from Ea. (|18f) that in tight-binding limit the in- 
stability opens at fc = Trh/2d similar to what happens 
in Bose-Einstein condensates Q • One should notice that 
results (|17I18(I , being derived from hydrodynamic theory, 
holds only for small values of q ■ There are however no 
restrictions on the value of the stationary velocity k/m 
of the fluid with respect to the lattice. 

For cf <C 2(5, where m = m* , the Doppler term is ex- 
actly the same as for Bose condensates. For higher values 
of the density, however, this term becomes smaller and 
vanishes identically for ep > 25, where to is linear in the 
density n. In this limit the current becomes density inde- 
pendent while the chemical potential no longer depends 
on the velocity of the fiuid. 

So far we have discussed the behaviour of a superfluid 
gas confined by a periodic potential. In the presence of 
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the additional harmonic potential Vho the low energy os- 
cillations exhibit new features. Of particular importance 
is the center of mass oscillation along the direction of the 
optical lattice. In the absence of the periodic potential, 
this mode would oscillate exactly at the frequency of the 
harmonic trap. The presence of the periodic potential 
give s rise to new interesting features. In a recent work 
|l2l | it has been shown that in a non-interacting Fermi gas 
with Fermi energy ep larger than the Bloch bandwith 26, 
the center of mass cannot oscillate around the bottom of 
the trap, but remains localized at one side of the har- 
monic field. Inclusion of collisions in a two-component 
gas favours the relaxation towards the equilibrium con- 
figuration. However, collisions cannot restore the propa- 
gation of sound in the hydrodynamic regime as happens 
in the absence of the periodic potential. In fact, it was 
concluded in Ref.0] that, under the condition ep > 2(5, 
the center of mass oscillations of a non superfluid gas are 
overdamped in the coUisional regime, due to umklapp 
processes. 

As recently suggested by Wouters et al. 14], in the 
superfluid regime the center of mass can instead oscillate 
and the corresponding frequency can be calculated start- 
ing from the hydrodynamic equations © and We 
will make the natural ansatz 



(j) — a{t)z, n — no{z — a{t),r±), 



(19) 



for the phase and for the density distribution where a 
and a are functions of time and no{r) — 2 J Q{ep{r) — 
e(k))(ik/(27r;i)3 is the equilibr ium density evaluated in 
the local density approximation with eF{r) — ej? — V/io(i")- 
This ansatz corresponds to a rigid shift of the density 
in coordinate space accompained by a uniform velocity 
field Vz = ha /2m. By inserting the ansatz H19|) into 
the Lagrangian (j2Jl and retaining only the lowest order 
(quadratic) terms in the functions a and a, one can recast 
in the form L(a, a, a) cx rabjia? 12+0^ l^mr;M+aaj2, 
where we have introduced the averaged effective mass |l5l | 



— = [ 4no(r)dr 
mcAi N J m 



N 



(de/dKf5{eF[v) - e{\i)) dkdr / {2tt hf (20) 



and N = J no{r)dr is the total number of particles. The 
above Lagrangian describes a simple harmonic oscillator 
with frequency 



mcM 



(21) 



Result (|21|l for the frequency of the center of mass os- 
cillation can be also derived microscopically using a sum 
rule approach, based on the ratio 



(M^ 



mi 
m_i 



(22) 




FIG. 1: Function / [See Ea. l|28|l ] versus the parameter y = 
eF/2S (solid line). The asymptotic limit 5/8y is also shown 
(dashed line). 



between the energy weighted and inverse energy weighted 
moments relative to the dipole operator D = J2j=i ^j- 
The energy weighted moment can be calculated using the 
general expression mi = {1/2){[D,[H,D]]) where () is 
the average over the ground state and the effective Hamil- 
tonianT? = J2j <^{^3) + J2j Vho{rj)+J2j<e^'^-body{\ r/c- 

I) is the projection of the full Hamiltonian into the low- 
est Bloch band. One finds mi = (l/2)('^^^^ d^e/dk^), 
which is the analog of the f-sum rule 16]. On the other 
hand the inverse energy weighted moment is easily re- 
lated to the dipole static polarizability api, fixed by the 
harmonic potential: m_i = (l/2)az3 = N/{2mul). By 
evaluating the average (X]^i ^^^/^^zj) the limit of a 
weakly interacting gas and inserting the results for mi 
and m_i into Ea. (|22|l . one recovers result H21|) for the 
frequency of the center of mass oscillation. 

Equations (|20|l . H21|) hold for any value of the laser 
intensity provided the system remains superfluid. In the 
tight binding limit the formalism becomes particularly 
simple. In fact in this case the effective mass (jJOJ, can 
be conveniently written as 1 /mc m — f I'm* , where 



b j h{xfl^ sm^ xQ[h{x)]dx 
■' ^ 2 f h{x)^/^e[h{x)]dx 



(23) 



is a dimensionless function of the ratio eF/26 and the de- 
pendence of m* = h?J6d^ on the laser intensity s is given, 
for example, in Ref.fg]- The function h{x) is defined by 
h{x) — ei?/(5— l+cosx. Here ep is the Fermi energy which 
is related to the free value = {3N)^^^h{LUA_^ujzV^^ 
evaluated in the absence of the periodic potential, by 
the equation 



0.3 ^}^(ErY/^ 



(4) 



57T^ \ S J 



h{xf/'^e[h{x)]dx, (24) 



which follows from the normalization condition N 
J no(r)dr. 
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In Fig^ we plot / as a function of the parameter y = 
ep/'^S. In the hmit ep 'ti 25, m does not depend on the 
density and one finds f — 1 and tticm = ™* — fi^/ScP. 
For higher values of the density, the frequency of the os- 
cillation becomes smaller because mcM increases with 
the density. In the limit ep ^ 25, we find the asymp- 
totic behavior / ~ bS/Aep (see dashed line in Fig^ and 
mcM = 'ih?ep/5(fS'^. 

As an example, we consider a two-component gas of 
N — 10^ potassium {'^'^K) atoms with trap frequencies 
uj± = 2tt ■ 275Hz and ujz = 2tt ■ 24Hz, corresponding 
to e^p/kB — SOOnif. For the optical lattice we assume 
s = 5 and periodicity d — 400nm, corresponding to Epi = 
9.25 = 37AnK. From Eq.(|23 one finds ep = 0.85e"p and 
from the value ep/2S = 4.1, one finds /(4.1) = 0.17, 
yielding mcM/'m = 10.94. Ea. (|21|l predicts a significant 
reduction {ujcm — O.SOuJz) of the frequency of the center 
of mass oscillation. Notice indeed that this reduction is 
larger than in the corresponding case of an oscillating 
Bose-Einstein condensate, where ujcm = y^m/m*uJz = 
0.73uz. 

In conclusion we have investigated the dynamic be- 
haviour of a superfluid Fermi gas in the presence of a ID 
optical lattice. We have calculated the velocity of sound 
and studied the conditions of stability (both energetical 
and dynamical) against the creation of elementary exci- 
tations when the fluid moves with respect to the lattice. 
Moreover, we have calculated the frequency of the center 
of mass oscillation (in the BCS limit) when the gas is 
further confined by a harmonic trap. Our results might 
be useful for the identification of the superfluid phase 
of interacting Fermi gases. Near a Feshbach resonance, 
where the scattering length becomes very large, one still 
expects the sytem to remain superfluid for moderate in- 
tensities of the laser field and, consequently, to exhibit 
undamped center of mass oscillations in the presence of 
harmonic trapping. In this case, however, the evaluation 
of the effective mass cannot be based on Eas. (|10l20|l and 
requires a more complete many-body calculation. 
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